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We perform population synthesis calculations of present-day post common 
envelope binaries (PCEBs) and zero-age cataclysmic variables (ZACVs) using 
a common envelope efficiency parameter, cxce, that is a function of secondary 
CN [ mass, M s . We investigate three basic possibilities: (1) a standard constant o>ce 
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prescription, with cxce = 1-0, 0.6, 0.3, 0.2, 0.1 and 0.05, to provide a baseline for 



comparison, (2) a power law dependence, acE = (M s ) n , with n = 0.5, 1.0 and 
\q [ 2.0, and (3) a dependence in which acE approaches 1 for large secondary masses 

and acE = below some assumed cutoff mass, cxce = 1 — M cut /M s , where M cut is 
the cutoff mass and is equal to 0.0375, 0.075 and 0.15 M . For each population, 



we compute orbital period, orbital separation, white dwarf mass and secondary 
mass distributions. 



We find that if acE ^ 0.2 in our constant cxce sequence, the predicted 
present-day ZACV population is significantly modified compared with our stan- 
dard model (cxce = 1-0). All prior population synthesis calculations of the for- 
mation of CVs only considered values of acE > 0.3 and found that their model 
populations were not strongly dependent upon the value of oice- Our results 
indicate that a much wider range of values for oice-, including very low values, 
must be considered in order for a dependence to be seen. In our variable acE se- 
quences for ZACVs, we find that for models in which cxce decreases very rapidly 
for small secondary masses, the orbital period distribution below the period gap 
differs significantly from our standard model. These differences are most evident 
in our power law sequence model with n = 2 and in our cutoff mass sequence 
model with M cut = 0.15 M Q . In these two models, the fraction of CVs forming 
with orbital periods below the gap is reduced significantly, the fraction form- 
ing in the gap is increased significantly, and both the short-period peak and the 
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minimum period in the ZACV orbital period distribution shift to considerably 
longer orbital periods compared with our standard model. We suggest that the 
observed scarcity of CVs with P < 77 min may possibly provide evidence that 
progenitor binaries with very low mass secondaries (M s < 0.10 M ) are unable 
to avoid merger within the common envelope. We also suggest that if acE de- 
creases rapidly for small secondary masses, as in our power law sequence with 
n > 1, it is possible that the lower edge of the period gap could be, in part, 
an imprint of the ZACV population. Such an imprint could be important in 
recently-proposed non-standard scenarios of CV secular evolution, such as cir- 
cumbinary disks, which have difficulty in reproducing the sharpness of the lower 
edge of the gap. 

In our constant cice model sequence for present-day PCEBs, we find that for 
all values of olce, the majority of the systems contain secondaries with masses > 
0.375 M , orbital periods > 1 day and orbital separations > 0.025 AU, with most 
having periods of ~ 3 days and separations of ~ 0.05 AU. These models further 
predict that the present-day population of PCEBs should contain roughly an 
equal number of systems with He and CO WDs if acE is globally near unity, but 
should be clearly dominated by systems containing CO WDs if acE is globally 
small (< 0.30). In our model sequences of present-day PCEBs in which acE is a 
function of secondary mass, the only distribution that varies significantly is the 
secondary mass distribution. In the power law model sequence, as n is increased 
from 0.5 to 2.0, the fraction of PCEBs with M s < 0.375 M decreases by a factor 
of ~ 4 from 0.26 to 0.06. In the n = 2.0 model, there are no present-day PCEBs 
with secondary masses less than 0.10 M Q . In the cutoff mass model sequence, 
significant changes only occur near the cutoff mass and the distributions are 
nearly identical for M s > 0.375 M . 

Based on the results of this investigation, we suggest to theorists who perform 
detailed hydrodynamical calculations of common envelope evolution that a se- 
quence of models with a fixed mass giant and very low mass secondaries, ranging 
from 0.3 M down to 0.05 M , should be given some priority. We further suggest 
to observers that a well-defined, statistically complete sample of PCEBs, partic- 
ularly with regard to very low mass secondaries, is crucially needed to provide 
tests of detailed models of common envelope evolution. 

Subject headings: binaries: close — novae, cataclysmic variables — stars: evolution- 
stars: low-mass, brown dwarfs — subdwarfs 
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Introduction 



Common envelope (CE) evolution is believed to play an essential role in the evolution 
of many close binary systems. In particular, forming short period (~ days) close binary 
systems from wide progenitor binaries (~ years) necessitates some phase of evolution in 
which a significant amount of the angular momentum is lost from the system. In a typical 
CE scenario, the primary star contacts its Roche lobe during a configuration in which mass 
transfer to the secondary star is dynamically unstable, such as when the primary is a giant 
and has a deep convective envelope. The secondary star cannot accommodate the accretion 
of the material at such high mass transfer rates , and instead becomes qui ckly engulfed in the 
envelope of the primary (see llben fc Liviolll993l ; lTaam fc Sandquistll2000l for recent reviews). 
The secondary and the core of the primary star are now orbiting within the envelope of the 
primary, which has become "common" to both components and which is not co-rotating 
with respect to the orbit of these components. Dynamical friction between the orbiting 
components and the non-corotating envelope cause a dissipation of the orbital energy of 
the components into the CE. If enough energy is able to be transferred from the orbit of 
the engulfed components to the envelope before the components merge, then the envelope 
may be able to become unbound and be shed. The result is a much shorter period binary 
consisting of the core of the primary and the original seco ndary, which exits the CE relatively 
unscathed (e.g.. IWebbinklll988l ; iHjellming fc Taamlll99ll ). 



Modeling of the CE phase has been undertaken by several authors and dramatic improve- 

ments in sophistication have been made over the past thirty years (e.g-. lMeyer fc Meyer- Hofmeister 
19791: a series of 11 papers by R. Taam, P . Bodenheimer and co-workers, the latest being 
Sandquist et al.ll2000l ; iRasio &: Livid 1 19961 ) . The most recent papers on thi s topic involve 



3-dimensional hydrodynamic calculations of this phase of evolution (e.g., IRasio &: Livio 



19961 : ISandquist et~alfll998l . I2OO0I : be Marco et~alTl2003al lbh . Despite these advances, the 
latest model calculations are not yet able to follow the evolution through full envelope ejec- 
tion and thus can on l y pla ce upper limits on the final (post-CE) orbital separations (see 
Sandquist et al.l Il998l . l2000l ) . This inability is primarily due to the large range in spatial 
scales that need to be resolved in following the evolution from start (orbital separation of 
~ 100-1000 R ) to finish (orbital separation of ~ a few R©). Consequently, the results of 
these detailed model calculations of the CE phase have not been able to be included into 
population synthesis calculations as yet. 

Instead, population synthesis calculations hav e relied on simple energy co nsiderations 
to relate the post- and pre-CE orbital parameters (ITutukov fc Yunglesonlll979l ). Although 



the exact prescription varies somewhat from author t o author, a typical ex 



pression used in 



population synthesis calculations is given below (e.g., Ilben fc Tutukovlll985l ; |Politanolll988 
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19961 : |Di Stefano fe Rappaportlll994l : IWillems et al J 12005( 1: 
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where M p is the total primary mass at the onset of the CE phase, M s is the mass of the 
secondary, M c is the core mass of the primary, M e is the envelope mass of the primary (M e 
= M p — M c ), Ri is the radius of the primary at the onset of the CE phase (and is therefore 
equal to the radius of the primary's Roche lobe), Oj and a/ are the pre-and post-CE orbital 
separations, respectively, A is a dimensionless factor related to the structure of the giant 
star, and occe is a parameter that measures the efficiency at which orbital energy is able to 
be used in unbinding the CE. Essentially, this prescription approximates the CE phase by 
saying that when some fraction of the orbital energy released in the spiral-in process equals 
the initial gravitational binding energjj^of the primary's envelope, the CE becomes unbound 
and is presumably (somehow) ejected. The parameter, ctcE, embodies a major uncertainty 
in this simplified prescription. 



Standard formation models of CVs have typically chosen otcR = 1.0 and/or pick = 0-3 



Politano 



2001 



1988 



Willems et al 



1996; de Koollll992l : iHan. Podsiadlowski fe Eggletonlll995l : iHowell et al. 



20051 ) . The latter choice reflects the few estimates of ccce from detailed 



models of the CE phase. These estimates suggest that acE is les s than 1.0 and may be closer 



to 0.3 to 0.5 , depending on the component masses of the binary (jTaamlll996l ; ISandquist et al. 



19981 . |2000| ; also see Fig. 1). Regardless of the exact value chosen for it, virtually all 



population synthesis calculations have assumed that olce is a global constant for all binaries 
in the population. However, it is unlikely that this is indeed the case. A more realistic 
assumption is that the value of olc e depends on t he orbital parameters of the binary, such 
as the mass of the secondary (see iPolitanol 120041 ) and the state of evolution of the giant 
(e.g., its core mass and total mass). In this series of papers, therefore, we undertake the 
first systematic investigation of the effect of a variable olce (i-e., one that depends on the 
orbital parameters of the binary) on population synthesis of close binaries. Specifically, in 
this first paper, we compute population models of present-day post-CE binaries (PCEBs) 
and zero-age cataclysmic variables (ZACVs) using a CE efficiency parameter that depends 
upon the mass of the secondary star. In future papers in this series, we will investigate 
separately the dependence of oice on other orbital parameters. 



1 We note that some authors include the thermodynamic internal energy of the envelope on the right side 
in eq. 1 (e.g., Han. Podsiadlowski fe Eggleton 19951 : Podsiadlowski. Han fe Rappaport 2003 ). The right side 
of eq. 1 would then represent the total binding energy of the envelope (gravitational potential energy + 
internal energy). Nevertheless, the meaning of acE as a parameter that measures the efficiency of energy 
transfer from the orbit to the envelope rema ins the same regardless of whether th e total binding energy or 
the gravitational binding energy is used fsee lHan. Podsiadlowski fe Eggletonlll995T ). 
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A central goal in undertaking the present study is to provide guidance to researchers 
who calculate detailed hydrodynamical models of CE evolution. Such models are highly 
computationally intensive (Taam, priv. comm. 2005). Population synthesis calculations, 
on the other hand, are relatively quick and are ideally suited for identifying interesting 
regimes within a specified parameter space that bear further, detailed investigation. In 
addition, identification of these regimes may suggest future observations to provide tests of 
the hydrodynamical models. 



Further motivation for the present investigation comes from recent work by iPolitano 



( 120041 ). who investigated the formation of cataclysmic variables (CVs) with brown dwarf 
(BD) secondaries. For standard assumptions (i.e., o>ce = 1 and a flat mass ratio distribution 
in zero-age main sequence [ZAMS] binaries), his models predict that 18% of the present-day 
population of ZACVs contain substellar secondaries. Unlike period-bounce CVs, where the 
secondary became substellar as a result of mass transfer to the WD, the secondaries in 
Politano's study were born substellar. What is particularly interesting about ZACVs with 
subs tellar secondar ies is that the majority of them (80%) have orbital periods shorter than 77 
mm jPolitanokool . the observed orbital period minimum in CVs. From his models, Politano 
predicts that as much as 15% of the present-day CV population could contain secondaries 
that were born substellar a nd have orbital perio ds less than 77 min. Yet of the ~ 600 CVs 
with known orbital periods (IRitter fc Kolbll2003l ) . only three have orbital periods less than 77 
min: V485 Cen (59 min), EI Psc (64 min) and J1507+5230 (67 min ). Moreover, it is highly 



unlikely that the secon d aries in V485 Cen an d EI Psc are B Ds (lAugusteijn et al.l 11996 



Thorstensen et al.l 120021 ; lUemura et al.ll2002l ). iPolitanol (120041 ) suggests that one possible 



reason for this large discrepancy between his models and observations is that progenitor 
binaries that contain BD secondaries may not be able to survive the CE phase. This would 
imply that olce is a function of secondary mass and that there exists a cutoff secondary mass 
below which merger is inescapable. This is not unreasonable, since the rate of dissipation of 
orbital energy within the CE will decrease as the secondary mass decreases. Consequently, 
below some secondary mass, this dissipation will be insufficient to drive hydrodynamical 
motion but rather will be carried away by convection within the CE (R. Taam, 2003, private 
communication). We examine the implications of Politano's suggestion in this paper by 
considering a functional dependence of olce on M s that contains a low-mass cutoff. 

In the next section, we describe the method used in our population synthesis calculations 
and the functional relationships for (xce(M s ) that we investigated. In section 3, we present 
our model results for the present-day populations of PCEBs and ZACVs. In section 4, 
we discuss our results and suggest future model calculations of the CE phase and future 
observations of PCEBs. We close with a summary of our main conclusions. 
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Method 



A Monte Carlo pop ulation synthesis code has been used in this study. This co de is the 
same as the one u sed in iPolitano &: Weilerl (120061 ) and the physics is described in iPolitano 
(119881 . Il996l . 120041 ) except as noted below. We refer the reader to those papers for detailed 
discussions, and here summarize only the key assumptions and features of the code and the 
modified treatment of the CE phase. 

In our Monte Carlo calculations, we begin with 10 7 ZAMS binaries. Following a standard 
approach, we assume that the distribution of ZAMS binaries is the prod uct of separate 
distri butions over primary mass, mass ratio and orbital period. We use a Miller fc Scalo 
Jl979h IMF for the primary masses, which we r eproduce numerically using the following 
Monte Carlo representation (jEggleton et al.lll989n : 



M p (u) 



0.19u[(l 



u 



.3/4 



+ 0.032(1 



u 



1/41-1 



(2) 



where u is a random number uniformly distribu ted betwe en and 1. We assume a ZAMS 
orbital period distribution that is flat in log P (Abt 1983 1^1 and a ZAMS mass ratio distri- 
bution that is flat in q (i.e., q(q) dq = I dq), where q = M s /M p (jDuquennoy &: Mayor! Il991 



Mazeh et al.lll992l ; iGoldberg et al.l 120031 ) . Once the primary mass and mass ratio have been 
generated for a given binary, the secondary mass is determined by their product: M s = q M p . 
The maximum primary mass in our calculations is 10 M . The minimum secondary mass 
that we consider is 0.013 M Q . 

To form a ZACV in our population synthesis code, a given ZAMS binary is evolved to 
the point when the primary contacts its Roche lobe while it is ascending the giant branc h (es). 
Wind loss during ascent of the giant branch(es) is incorporated via a iReimers! (119751 ) pre- 
scription. Relationships between the radius of the giant, its core mas s and its total mass 
are derived from analytic fits to detailed stellar evolution models (see iPolitand Il988l . 11996 
for an in-depth discussion and for detailed references). Simple energetic considerations are 
still used to relate the pre- and p ost-CE orbital separ ations for the CE phase (see eq. 1). 
Consistent with other studies (e.g.. IWillems et al.ll2005l ). we choose A = 1. However, we now 
allow acE to be a function of secondary mass (see below). Since the CE phase is so brief 



We no te that such a distribution i s inconsistent with the orbital period distributions found by I Griffin 



() 19851 ) and lDuquennov fc Mavorl (|l99lf ). However, the use of a ZAMS period distribution tha t is flat in log 



P (or an orbital separation distribution that is flat in log a) is widespread in the literature (e.g.. lHurlev et al 



2002t IWillems et al.l 12005c IZhang et al.l 120051 ) and thus facilitates comparison. Additionally, model pop- 
ulations of ZACVs and PCEBs that we calculated using a Griffin or Duquennoy and Mayor-type period 
distribution produced substantially the same conclusions as those calculated using a period distribution that 
is flat in log P. 
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(e.g., Ilben fc Livid Il993l ; iTaam fc Sandquistll2000l ). the mass of the core remains essentially 
unchanged during the CE phase and, thus, the white dwarf (WD) mass upon exit of the CE 
phase is equal to the core mass of the giant at the onset of the CE ph ase. Lastly, the clas- 
sical prescription for gravitational radiation ( Landau &: Lifshitzl 1951) and the prescription 
for disrupted magnetic braking (7 = 2) from lRappaport. Verbunt fc Joss! (119831 ) are used to 
describe angular moment um losses during the post-CE phase. We use detai l ed stellar mod- 



els fr om the Lyon group (IChabrier fc Baraffdll997l ; iBaraffe et al.Hl 9981 . 120031 ; IChabrier et al. 
2 OOOl ) for secondaries with masses less than ~ 0.5 M Q , including substellar secondaries. For 
secondaries with masses greater than ~ 0.5 M Q , we use fits to stellar models from Webbink 



Politanolll988l . Il996l ). In our models, the secondary is fully convective for M s < 0.37 



(see 
M . 

In addition to ZACVs, we model the present-day population of PCEBs. To avoid 
ambiguity, we note that when we use the term "PCEBs" in this paper, we are specifically 
referring to binaries that (1) are detached, (2) contain a WD primary and a MS or BD 
secondary, and (3) have undergone a single CE phase. Clearly, other important types of 
close binaries exist that have undergone a single CE phase-for exam ple, systems containing 



a giant secondar y and a WD, such as symbiotic stars or Ba stars (e.g.. lYungleson et al.lll995 



Han et al.lll995l ). However, we have limited the scope of our calculations to only PCEBs that 
satisfy the above three criteria. 

The same code as described above for ZACVs was used for PCEBs, but with two 
important modifications: 1) the secondary radius must be smaller than its Roche lobe radius 
at the present epoch and 2) we do not impose an upper limit on the binary's mas s ratio in 
order to satisfy stability against thermal or dynamical timescale mass transfer (see iPolitano 
19961 ). since the system remains detached at the present epoch. In calculating the population 
of present-day PCEBs, four timescales are of interest: (1) tj,, the time that the progenitor 
binary was formed (measured from the beginning of the Galaxy), (2) t eVtP , the time it takes 
the primary to evolve off of the main sequence, become a giant, and contact its Roche lobe 
to initiate the CE phase, (3) tpcEB, the time from the end of the CE phase until the present 
epoch, and (4), tc a u the age of the Galaxy. These four timescales must satisfy the following 
constraint, t b + t ev>p + t PC EB = ^Gai- We make the following assumptions regarding these 
time scales: (1) the stellar birth rate throughout the Galaxy's history has been constant, (2) 
the CE phase is so rapid t hat the time spent in it by the binary is neg ligible compared to 
the other time scales (e.g., Ilben fc Livid 1 19931 : ITaam fc Sandquistll2000l ). and (3) the age of 
the Galaxy is 10 10 yrs. 



For both populations, we have now allowed acE to be a function of secondary mass. We 
do not change the formal structure of eq. 1, but simply replace ace by a specified function of 
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secondary mass. There are only a handful of values of olc f. that have been dete r mined from 



detailed hydrodynamical calculations of the CE phase (jSandquist et al.l Il998l . l2000l ). and 
even these are only upper limits since Sandquist et al. were not able to follow the evolution 
completely through envelope ejection because of spatial resolution problems. Consequently, 
detailed models of the CE phase are able to provide little guidance in constructing an assumed 
dependence of q>ce on M s . We therefore chose two simple functional forms: (1) a power law 
dependence, oice = (M s ) n , with n = 0.5, 1.0, and 2.0 and M s in solar masses, and (2) a 
dependence in which olce approaches 1 for large secondary masses and is zero below some 
assumed cutoff mass, o>ce = 1 — M cut /M s , where M cut is the cutoff mass. We set M cut = 
0.0375, 0.075, and 0.15 M . These cutoff masses represent 1/2, 1, and 2 times the substellar 
transition mass, repectively. Figure 1 shows the acE estimates from detailed hydrodynamical 
models (x's in the figure) along with our two assumed functional relationships for o:ce(M s ). 
As seen in the figure, our assumed relationships span the parameter space containing the 
estimates of otcE from the detailed models. 

In addition to the variable olce models, we also calculated constant olce models for 
six values of olce'- 1-0, 0.6, 0.3, 0.2, 0.1, and 0.05. These constant acE models provide 
a baseline for comparison. In this paper, we shall refer to the olce — 1-0 model as our 
"standard" model. 

To avoid adding further uncertainty to our investigation, we have not normalized our 
model distributions to provide absolute numbers of PCEBs or ZACVs. Instead, our results 
will be described comparatively and will typically fall into one of two categories: a comparison 
of relative features within a given model distribution or a comparison of a particular feature 
across a given model sequence or across all model sequences. 



3. Results 

In this section, we present our model populations for present-day PCEBs and ZACVs 
calculated using the o:ce(M s ) relationships described in the previous section. For each model 
population, four distributions were calculated: orbital period, orbital separation, secondary 
mass, and WD mass. Our model populations of present-day PCEBs are discussed in section 
3.1 and our model populations of present-day ZACVs are discussed in section 3.2. Selected 
results for each model are given in Table 2 for PCEBs and in Table 3 for ZACVs. In these 
tables, the constant ccce models are designated by "CON" and the value of acE, the power 
law models are designated by "PL" and the value of the exponent n, and the cutoff mass 
models are designated by "CUTM" and the value of M cut . 
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Before proceeding to the specific populations, we discuss the fraction of systems that 
survive the CE phase, f surv , for each model. We define this fraction as the number of systems 
that escape merger during the CE phase divided by the total number of systems that entered 
the CE phase. We note that for our initial population of 10 7 binaries, ~ 17% of them enter a 
CE phase at some point during the Galaxy's history due to the primary filling its Roche lobe 
as a giant. Values of f surv are listed in Table 1. For the CON model sequence, we find that 
fsurv decreases roughly as (oice) ' 5 from 0.50 to 0.07 as oice is decreased from 1.0 to 0.05. 
Interestingly, the fraction of systems that survive CE evolution does not vary significantly 
in either the PL or the CUTM model sequences. In fact, if the CON 1.0 model is included 
as n = in the PL sequence, then over a wide range of dependences from olce independent 
of M s (n = 0) to o>ce very strongly dependent on M s (n — 2), we find that f surv decreases 
only slightly from 0.50 to 0.39. For the CUTM model sequence, a similar trend is found, 
with f surv decreasing from 0.48 to 0.43 as M cut is increased from 0.0375 to 0.15 M . 



3.1. PCEBs 

3.1.1. o>ce = constant 

The distributions of the orbital periods (top panels) and the orbital separations (bottom 
panels) in present-day PCEBs for acE = 1-0, 0.6, and 0.3 (left panels) and for ccce = 0.2, 
0.1, and 0.05 (right panels) are shown in Figure 2. The distributions of the secondary masses 
(top panels) and the WD masses (bottom panels) in present-day PCEBs for the same choices 
of olce are shown in Figure 3. 

Qualitatively, we find that decreasing cxce from 1.0 to 0.05 yields little change in the 
overall shapes of the orb ital period, orbital separa tion and secondary mass distributions for 



PCEBs. As discussed in lPolitano fc Weilerl (120061 ). the sharp drop in the number of present 



day PCEBs with secondary masses greater than 0.37 M Q is due to the increased efficiency of 
magnetic braking compared with gravitational radiation in shrinking the orbit and bringing 
systems into contact by the present epoch. Consequently, only PCEBs that exit the CE 
phase with relatively wide orbits (i.e., a > 0.025 AU or P > 1 day) remain detached at the 
present epoch. This is reflected in the large increase in the number of PCEBs near P = 1 
day and a = 0.025 AU for all values of acE in Fig. 2. 

For the purpose of comparison and to expedite the discussion, throughout this paper 
we will denote PCEB secondaries with M s < 0.375 M Q as "low-mass" secondaries and those 
with M s > 0.375 M as "high-mass" secondaries (the choice of 0.375 rather than 0.37 has to 
due with how the distribution was binned). Similarly, throughout this paper we will denote 
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present-day PCEBs with P < 1 day or a < 0.025 AU as "short-period" or "close" PCEBs 
and those with P > 1 day or a > 0.025 AU as "long-period" or "wide" PCEBs. 

Quantitatively, we find that as qzce is decreased from 1.0 to 0.05, the fraction of present- 
day PCEBs with low-mass secondaries decreases slightly from ~ 1/3 to 1/5. The fraction 
of present-day short-period PCEBs varies slightly, increasing from 0.23 for acE = 1-0 to a 
maximum of 0.27 and then decreasing to 0.19 for acE = 0.05. For all values of acE, the 
orbital period distribution peaks at ~ 3 days and the orbital separation distribution peaks 
at ~ 0.05 AU. Thus, over a wide range of values for o>ce in our CON sequence, our models 
predict that the majority of present-day PCEBs should contain high-mass secondaries, have 
long orbital periods and wide orbits, with most having periods of ~ 3 days and separations 
of ~ 0.05 AU. 

The distribution that is most affected by decreasing a CE is the WD mass distribution. 
As a C E is decreased from 1.0 to 0.05, the fraction of He WDs {M WD < 0.5 M ) in PCEBs 
decreases by a factor of ~ 20, from 0.51 to 0.02, and the fraction of CO WDg§ (M WD > 0.5 
M ) increases correspondingly from 0.49 to 0.98. Therefore, our CON model sequence pre- 
dicts that the present-day population of PCEBs should contain roughly an equal number of 
systems with He and CO WDs if q>ce is globally near unity, but should be clearly dominated 
by systems containing CO WDs if o>ce is globally small (< 0.30). 

Lastly, we find that the predicted relative number of present-day PCEBs compared 
with our standard model decreases monotonically as acE is decreased in our CON model 
sequence. For cuce — 0.05, the number of PCEBs is approximately an order of magnitude 
smaller than in our standard model. 



3.1.2. a CE = (M s ) n 

The distributions of the orbital periods (top left), the orbital separations (top right), 
the secondary masses (bottom left) and the WD masses (bottom right) in PCEBs at the 
present epoch for olce = {M s ) n , where n = 0.5, 1.0, and 2.0, are shown in Figure 4. 

Qualitatively, we find little change in the overall shapes of the distributions as n is 
varied, except for the secondary mass distribution. Quantitatively, the dominance of long- 
period systems is even stronger in the PL model sequence than in the CON model sequence. 
The fraction of short-period PCEBs is 0.19 for n = 0.5 and decreases to 0.08 for n = 2.0. 



3 Our code does not separately calculate the formation of PCEBs or ZACVs with ONeMg WDs. Instead, 
we assume that some fraction of CVs with high mass WDs contain ONeMg WDs. 
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As in the CON models, the distributions over orbital period and separation peak at ~ 3 
days and ~ 0.05 AU, respectively. The WD mass distribution is not affected as strongly 
in the PL models compared with the CON models. As n is increased from 0.5 to 2.0, 
the fraction of He WDs decreases slightly from 0.45 to 0.34 and the fraction of CO WDs 
correspondingly increases from 0.55 to 0.66. Thus, in our PL model sequence of present-day 
PCEBs, the majority of the systems have long orbital periods, wide orbits and contain CO 
WDs regardless of whether the dependence of acE on M s is weak or strong. 

The secondary mass distribution is strongly affected by the choice of n. As n is increased 
from 0.5 to 2.0, the fraction of PCEBs with low-mass secondaries decreases by a factor of 
~ 4, from 0.26 to 0.06. The decrease in the number of low-mass secondaries is particularly 
severe for the n = 2.0 model. In this case, there are no present-day PCEBs with secondary 
masses less than 0.10 M . 

The predicted relative number of present-day PCEBs compared with our standard model 
is 0.80 for n = 0.5 and decreases to 0.50 for n = 2.0. 



3.1.3. a CE = 1 - M cut /M s 

The distributions of the orbital periods (top left), the orbital separations (top right), 
the secondary masses (bottom left) and the WD masses (bottom right) in PCEBs at the 
present epoch for cxce = 1 — M cut /M s , where M cut = 0.0375 M Q , 0.075 M , and 0.15 M , are 
shown in Figure 5. 

The orbital period, orbital separation, and WD mass distributions are not affected 
significantly as the cutoff mass is varied. Inspection of Table 2 shows that the fraction of 
systems with He WDs, the fraction with short orbital periods and the fraction with small 
orbital separations each change only slightly as the cutoff mass is increased from 0.0375 to 
0.15 M . Further, the values of these fractions are not very different from those in our 
standard model. 

The secondary mass distribution is affected, but the effect is much more localized than 
it was for the PL model sequence. Significant changes occur only for low-mass secondaries 
and the distributions are nearly identical above ~ 0.375 M . The fraction of present-day 
PCEBs with low-mass secondaries decreases from 0.33 to 0.21 as the cutoff mass is increased. 
For M cut = 0.0375 M , the fraction with substellar secondaries is 0.007. This is about a 
factor of 4 smaller than in our standard model. The fraction with substellar secondaries is 
obviously zero for M cut = 0.075 M or 0.15 M . 
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We find that the predicted relative number of present-day PCEBs compared with our 
standard model decreases modestly from 0.92 to 0.72 as the cutoff mass is increased from 
0.0375 to 0.15 M R . 



3.2. ZACVs 

3.2.1. o>ce = constant 

The distributions of the orbital periods (top panels) and the orbital separations (bottom 
panels) in ZACVs at the present epoch for o>ce = 1-0, 0.6, and 0.3 (left panels) and for o>ce 
= 0.2, 0.1, and 0.05 (right panels) are shown in Figure 6. The distributions of the secondary 
masses (top panels) and the WD masses (bottom panels) in ZACVs at the present epoch for 
the same choices of cice are shown in Figure 7. 

Varying olce from 1.0 to 0.3 (left panels in Figures 6 and 7) in our CON model sequence 
has qualitatively little effect on any of the distributions. The fraction of systems containing 
He WDs decreases slightly from 0.58 to 0.46 and there is a slight increase in the fraction of 
systems with orbital periods above 2.25 hr (in and above the period gap). These results are 
in agreement with previ ous population synthesis calculations of ZACVs that chose o>ce = 



1.0 and a CE = 0.3 (e.g. Jde Koolll992ft . 



However, all distributions show marked differences as cice is varied from 0.2 to 0.05 
(right panels in Figures 6 and 7). From Table 3, the predicted fraction of CVs formed: 
below the period minimum (P < 77 min) decreases by a factor of 3 from 0.15 to 0.05; below 
the period gap (P < 2.25 hr) decreases by a factor of 2 from 0.53 to 0.27; in the period gap 
(2.25 hr < P < 2.75 hr) increases from 0.11 to 0.16, and above the period gap (P > 2.75 hr) 
increases from 0.35 to 0.56. The predicted fraction of CVs formed with low-mass secondaries 
decreases from ~ 2/3 to 1/2. The predicted fraction of CVs forming with He WDs decreases 
by greater than a factor of 3 from 0.41 to 0.12 and, correspondingly, the fraction of CVs 
forming with CO WDs increases from 0.59 to 0.88. 

Therefore, in our CON model sequence of present-day ZACVs, if o>ce <i 0.30, the 
majority of the systems contain He WDs and have orbital periods below the period gap. On 
the other hand, if o>ce ^ 0.20, then the majority of the systems contain CO WDs and there 
is a significant shift in the orbital period distribution towards periods in and above the gap. 
The majority of ZACVs contain low-mass secondaries regardless of the value of cxce- 

Lastly, we find that the predicted number of present-day ZACVs compared with our 
standard model decreases monotonically as o>ce is decreased in our CON model sequence. 
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For acE = 0.05, the number of ZACVs is approximately a factor of 4 smaller than in our 
standard model. 

3.2.2. a CE = {M s ) n 

The distributions of the orbital periods (top left), the orbital separations (top right), 
the secondary masses (bottom left) and the WD masses (bottom right) in ZACVs at the 
present epoch for acE = (M s ) n , where n = 0.5, 1.0, and 2.0, are shown in Figure 8. 

Unlike PCEBs, all ZACV distributions are affected significantly by the choice of n, and 
these effects are not restricted solely to n — 2. In Table 3, we see that as n is increased 
from 0.5 to 2.0, the fraction of CVs formed with orbital periods: less than 77 min is reduced 
from 0.15 to zero; below the period gap decreases by a factor of ~ 7 from 0.58 to 0.09; in 
the period gap increases by almost a factor of 3 from 0.07 to 0.18; and above the period gap 
doubles from 0.35 to 0.72. The fraction of CVs formed with low-mass secondaries decreases 
by almost a factor of 2 from 0.67 to 0.35. For n = 2, no CVs are formed with secondary 
masses less than ~ 0.10 M Q . The fraction of CVs forming with He WDs decreases by almost 
a factor of 6 from 0.46 to 0.08 as n is increased from 0.5 to 2.0 and, correspondingly, the 
fraction of CVs forming with CO WDs increases from 0.54 to 0.92. 

Thus, in our PL model sequence of present-day ZACVs, the majority of the systems 
contain He WDs regardless of whether the dependence of acE on M s is weak (n = 0.5) or 
strong (n = 2.0). Systems with low-mass secondaries dominate if the dependence of acE 
on M s is weak to moderate (n < 1), whereas systems with high-mass secondaries dominate 
if the dependence is strong. Finally, systems with orbital periods below the period gap are 
favored if the dependence of acE on M s is weak, while systems in and above the period gap 
are favored if the dependence is moderate to strong. 

The relative number of ZACVs compared with our standard model is 0.84 for n = 0.5 
and decreases to 0.37 for n = 2.0. 

3.2.3. a CE = l-M cut /M s 

The distributions of the orbital periods (top left), the orbital separations (top right), 
the secondary masses (bottom left) and the WD masses (bottom right) in ZACVs at the 
present epoch for a C E = 1 — M cut /M s , where M cut = 0.0375, 0.075, and 0.15 M Q , are shown 
in Figure 9. 
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Again, unlike PCEBs, all ZACV distributions are affected by the choice of M cut , although 
the effects are significant only for the lower end of each distribution (i.e., short orbital 
periods/separations and low-mass secondaries/WDs). We find that as the cutoff mass is 
increased from 0.0375 M to 0.15 M in our CUTM sequence, the fraction of CVs formed 
with orbital periods: less than 77 min decreases from 0.18 to zero; below the period gap 
decreases by a factor of 2 from 0.62 to 0.31; in the period gap nearly triples from 0.05 to 
0.13; and above the period gap increases from 0.34 to 0.55. The fraction of CVs formed with 
low-mass secondaries decreases from 0.68 to 0.48. This decrease is similar to the decrease 
in the fraction of low-mass secondaries observed in the PL model sequence. The overall 
number of CVs forming with He WDs decreases, while the overall number of CVs forming 
with CO WDs remains approximately the same (see Fig. 9), causing the fraction of present- 
day ZACVs with He WDs to decrease by almost a factor of 2 from 0.53 to 0.28 and the 
fraction with CO WDs to increase correspondingly from 0.47 to 0.72. 

Thus, in our CUTM model sequence of present-day ZACVs, the majority of the systems 
contain low-mass secondaries and have orbital periods below the period gap if M cut < 0.075 
M , while the majority of the systems contain high- mass secondaries and have orbital periods 
above the gap if M cut = 0.15 M Q . Further, approximately the same number of systems 
contain He WDs and CO WDs if M cut < 0.075 M , while the majority of the systems 
contain CO WDs if M cut = 0.15 M . 

Because the secondary star fills its Roche lobe in a ZACV, unlike in a PCEB, there is 
a strong correlation between the mass of the secondary and the orbital period in a ZACV. 
Consequently, the existence of a cutoff mass for the secondary implies the existence of a 
corresponding minimum orbital period in the ZACV orbital period distribution. For M cut = 
0.0375 M , there are no systems with the orbital periods shorter than 43 min. This minimum 
period is the same as in our standard model. However, when M cut is increased to 0.075 M , 
the minimum period increases to 55 min (see Fig. 9). Finally, when M cut is increased to 
0.15 M , the cutoff period increases to 102 min, over twice the value in our standard model. 
These results are significant since they clearly imply that if a cutoff mass for the secondary 
indeed exists and is greater than the substellar transition mass, the standard model for the 
formation of CVs below the period gap will need to be modified (see section 4.2). 

Lastly, we find that the relative number of ZACVs compared with our standard model 
decreases from 0.86 to 0.52 as the cutoff mass is increased from 0.0375 to 0.15 M . 
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Discussion 



4.1. Observations of PCEBs: the current state 



4-1.1. Large Surveys 



The observed number of detached WD-MS systems continues to increase dramatically. 
As few as six years ago, there were only ~ 40 known (Hillwig et al. 2000). That number 
has risen to over 800, mainly due to searches for them in large observational surveys such 
as the Sl oan Digital Sky Surve y (SDSS) and the Two Micron All Sky Survey (2MASS). For 
example, ISilvestri et al.l (120061 ) have recently compiled a catalog of 746 spectroscopically- 
i dentified, detached close binary systems in the SDSS through the Fourth Data Release 
( lAdelman-McCarthy et al.ll2006l ). The vast majority of these systems (~ 700) are detached 
WD-M dwarf binaries. Radial velocity follow-up studies have been performed and orbital pe 
riods are being determined in ~ 25 systems (N. Silvestri 2 006, priv. comm.) . IWachter et al. 



(120031 ) have searched the 2MASS point source catalog (jCutri et al.l 120031 ) for WDs with 
infrared excess, which may indicate the presence of an unresolved red dwarf companion, 
and identified 95 systems. They have followed up these identifications with Hubble Space 
Telescope (HST) snapshot observations to determine if i ndeed these WDs have low mass 
companions. First results from these HST observations (IFarihi. Hoard, fc Wachterl 120061 ) 
indicate that of the 48 systems studied thus far, 27 have fully or partially resolved com- 
panions and 15 more are almost certainly unresolved binaries. The resolved systems have 
estimated projected orbital separations between ~ 10 - 1000 AU, making it unlikely that any 
of these systems underwent a CE phase. Esti mated upper limits for the unreso lved binaries 
suggest orbital separations less than ~ 5 AU. IFarihi. Hoard, fc Wachterl (120061 ) suggest that 
these unresolved systems are likely WD + red dwarf pairs in close orbits and may be good 
candidates for radial velocity variables. Assuming the 48 systems studied thus far are repre- 
sentative of the entire sample of 95, they expect to resolve ~ 55 wide systems and identify 
a total of ~ 30 unresolved close binaries for future radial velocity studies. 



4-1.2. Compilations 



Schreiber fc Gansickd (120031 ) compiled a sample of 30 well-studied PCEBs with known 



orbital periods. Since they were primarily interested in PCEBs that would become CVs 
within a Hubble time, they restricted their sample to systems that have mass ratios less 
than one (i.e., M s < M\y D ), have orbital periods shorter than 2 days, and have main se- 
quence (or BD) secondaries (i.e., no subgiant or giant secondary stars). More recently, 
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Morales-Rueda et al.l (120051 ) compiled a list of all detached binaries containing at least one 
WD and having known orbital periods. In their compilation, 32 binaries are confirmed WD 
+ red dwarf pairs and 9 are sdOB + red dwarf pairs. Their sample includes 28 of the sys- 
tems in the ISchreiber fc Gansickd (120031 ) compilation. To the best of our knowledge, the 
Morales-Rueda et al.l J2005h compilation provides the most extensive list of PCEBs with 
known orbital periods that is currently available, although once orbital periods are de- 
termined for the subset of systems from SDSS and 2MASS, that number should increase 
significantly. 



Regrettably, neither the PCEBs in the iMorales-Rueda et al.l (120051 ) compilation nor 
those culled from the large SDSS and 2MASS data releases provide a well-defined, sta- 
tistically complete s ample of PCEBs. For example, the majority of the systems in the 
Silvestri et all (I2006h SDSS sample were targeted for sp ectroscopy because th eir colors re- 
sembled other, higher priority objects such as quasars (see iRichards et al.ll2002l for additional 
details). In addition, PCEBs containing either a hot WD and a late M/L type secondary or 
a cool WD and an early M/late K type seco ndary present particular obs ervational challenges 
and have likely been missed generally (see ISchreiber fc Gansickd 120031 for a nice discussion 
of this). In the former case, the WD is much brighter than the secondary in the optical and 
the system may be mistaken for a single WD based on its colors. In the latter case, the 
WD is hidden by the glare of the secondary and the system may be misidentified as a single 
early M or K dwarf based on its colors. A proper sample of PCEBs for comparison with 
population models such as those presented in this paper does not exist as yet. Until it does, 
any comparison with theoretical populations of PCEBs would be premature. However, the 
rapid growth in the number of observed PCEBs over just the past few years suggests that 
such a sample may be feasible within the next several years. Efforts are already in progress 
to improve the observed sample of PCEBs over the next 2-3 years (B. Gansicke, priv. comm. 
2006). 



4.2. Observations of CVs 



In the latest online version (RKcat7.6, 1 January 2006) of the lRitter fc Kolbl (120031 ) cat- 
alog, there are 604 CVs with known orbital periods (we exclude systems that RK designate 
as possible AM CVn stars). The distribution of orbital periods in observed CVs provides, 
by far, the best test of CV population models. We show this distribution in Figure 10. It 
is important to remember that the CVs in Figure 10 were born throughout the Galaxy's 
history and therefore do not represent present-day ZACVs. Consequently, because of sec- 
ular evolution, the distribution of orbital periods shown in Fig. 10 is most likely different 
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than the distribution of orbital periods in ZACVs (see e.g., 



Kolb 


1993; 


Howell et al. 


2001) 



Nevertheless, there are some features of the observed orbital period distribution in CVs that 
may provide clues as to the ZACV orbital period distribution and, therefore, are germane 
to the present discussion. 



4-2.1. The scarcity of CVs with orbital periods shorter than 77 min: evidence for a cutoff 

mass? 



One of the most prominent features in the observed orbital period distribution of CVs 
is the existence of a minimum period at 77 min. Of the ~ 600 CVs in this distribution, 
only 3 have orbital periods shorter than 77 min: V485 Cen (59 min), EI Psc (64 min), 
and J1507 +5230 (67 min). According to the accepted theoretical explanation of the period 
minimum (jPaczyhski fc Sienkiewiczl Il98ll ; iRappaport. Joss fc Webbink! Il982h , CVs evolve 
to shorter orbital periods until at some point the secondary has lost so much mass that 
its thermal timescale becomes longer than the mass transfer timescale. This point roughly 
coincides with the secondary having lost so much mass that it has also become substellar. 
Once this point is reached, the secondary's response to mass loss changes and it now expands 
upon mass loss rather than con tracts. Further mass transfer then results in an increase in 
the o rbital period of the system (jPaczyhski fe Sienkiewiczlll98ll ; IRappaport. Joss fc Webbink 
19821 ). Consequently, CVs that form with orbital periods longer than the orbital period 
minimum cannot evolve to orbital periods shorter than it. 

This explanation of the period minimum is suitable for CVs that form with orbital pe- 
riods greater than the period minimum, b ut it does not ad dress CVs that form with orbital 
periods shorter than the period minimum. iPolitanol (120041 ) performed the first detailedpop- 
ulation s ynthesis c alcula tions of the formation of CVs with ultrashort orbital periodso As 
noted in IPolitanol (120041 ). the secondaries in CVs that form with orbital periods less than 77 
min are very low mass or substellar. Consequently, realistic population synthesis calculations 
of the formation of ultrashort-period CVs awaited the development of re liable detailed mod- 
els of very low mass stars and BDs over the past five to ten years (e.g., IChabrier fc Baraffe 
2000h . 



Our s tandard model u ses essentially the same input physics for very low mass/substellar 



stars as in 



Politand (120041 ) (see section 2), and it predicts that almost 1/4 of all present-day 



4 We note that AM CVn stars have ultrashort orbital periods, but these systems have helium-rich donors. 
By ultrashort-period CVs, we are specifically referring to systems containing hydrogen-rich donors that have 
orbital periods below 77 min. 
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ZACVs form with orbital periods shorter than 77 min (see Table 3 and Figure 6). Where 
are these systems? Or to rephrase the question, how do we explain the current lack of CVs 
with orbital periods shorter than 77 min? We suggest four possible answers: (1) all CVs ever 
formed below 77 min have secularly evolved to orbital periods longer than 77 min by the 
present epoch; (2) CVs formed below the period minimum are too faint to see in quiescence 
and have recurrence lifetimes too long to have seen an outburst in recent history; (3) the 
assumption in our standard model that oice = 1 is not valid; and (4) the assumption in our 
standard model that g(q) is flat is not valid. 



The first possibility is unlikely. iKolb fc Baraffd (119991 ) calculated secular evolution se- 
quences for CVs that formed with a 0.6 M WD and secondary masses between 0.04 and 0.21 
M Q using full stellar and substellar models for the secondary. These CVs formed at orbital 
periods between 54 min and 2.1 hr. For a 0.07 M secondary, they found that it took 1.5 
Gyr for the system to evolve from 54 min to 77 min@ Given that almost 1/4 of all CVs are 
predicted to form with orbital periods less than 77 min at the present epoch, it is difficult to 
understand how subsequent secular evolution could evacuate the entire orbital period regime 
shorter than 77 min. As a comparative example, a CV that formed at the upper edge of the 
period gap would t ake ~ 1 Gyr to cross the gap, assuming only gravitational radiation (e.g., 



Howell et al.lll997l ). Yet the period gap is not completely devoid of CVs. 



The second possibility is also unlikely. IKolb fc Baraffd (Il999l ) also computed mass 



transfer rates as a function of orbital period for the same evolutionary sequences discussed 
in the previous paragraph. They found that the evolutionary tracks merge into a common 
track at an orbital period of ~ 70 minutes regardless of whether the CV was formed with 
an or bital period longer than 77 min or shorter than 77 min (see Fig. 1 in IKolb fc Baraffe 



19991 ). Indeed, their models predict that present-day CVs with orbital periods less than 77 



min should have mass transfer rates greater than or equal to present-day CVs with orbital 
periods slightly greater than or equal to 77 min. Well known examples of the latter are 
WZ Sge and AL Com (P = 81 min for both systems). These two systems and other WZ 
Sge-type systems were not discovered in quiescence, but were discovered because of their 
outbursts. Consequently, if dwarf nova outbursts have been observed in these systems, then 
it is difficult to understand why similar outbursts would not have been observed in systems 
just below the period minimum, assuming they exist. Since such systems are predicted to 
have mass transfer rates similar to or greater than WZ Sge-type systems ( IKolb fc Baraffe 



19991 ). they should have recurrence times that are similar or even shorter. 



5 In the lKolb fc Baraffd (|l999l ) study, gravitational radiation was the only angular momentum loss mech- 
anism driving secul ar evolution. If an additional sourc e of angular momentum loss is allowed to operate 
below the gap (e.g., Pattersonl 1998 : Willems et al. 2005 ). this timescale would be reduced. 
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Determining the viability of the third possibility is one of our goals in undertaking the 
present investigation. To that end, inspection of Table 3 shows that the fraction of ZACVs 
with orbital periods shorter than 77 min is zero in 2 of the 12 models: PL 2.0 and CUTM 
0.15. A common feature of these models is that olce is exactly or essentially zero below 
some secondary mass. For the CUTM 0.15 model, that mass is 0.15 M . For the PL 2.0 
model, that mass is ~ 0.10 M Q . Consequently, these two models demonstrate that if merger 
within the CE is inescapable below a secondary mass of ~ 0.1 M Q or greater, agreement 
between the predicted and the observed orbital period distribution in CVs below 77 min 
can be greatly improved. We note that in 3 of the models, CON 0.05, PL 1.0 and CUTM 
0.075, the fraction of CVs formed below 77 minutes is 0.05. While this is still an order of 
magnitude higher than the observed fraction (at best, 3 out 604 or ~ 0.005), these may be 
viable models if subsequent secular evolution can move systems to periods longer than 77 
min on a sufficiently short time scale. 

The fourth possibility is that the distribution of ZAMS binaries containing a MS pri- 
mary and a very low mass/BD secondary is different than the distribution of ZAMS bi- 
naries containing two MS stars. Observations of companions to solar-type stars (spec- 
tral types F-M) indicate a relat ive scarcity of BP com panions compared with either less 
massive p lanetary companions (iMarcv fc Butler! 12000 or m ore massive stellar compan- 



ions (e.g., 



Duquennoy fc Mayor! Il99ll ; iFischer Sz Marcylll992l ). This comparative deficit of 
BP companions to solar-typ e stars has been termed the "brown dwarf desert." Recently, 



Grether &: Lineweaverl (120061 ) have described quantitatively the characteristics of the brown 
dwarf desert for nearby binaries (d < 50 pc) with orbital periods less than 5 yr (orbital 
separations < 3 AU). They find that ~ 16% of these nearby binaries contain solar type-stars 
with companions more massive than Jupiter: ~ 11 % h ave stella r comp anions, < 1% have 
BP companions, and ~ 5% have planetary companions. iPolitanol (120041 ) found that all ZA- 
CVs with BP secondaries evolved from progenitor ZAMS binaries with orbital separations 
less than 3 AU and ~ 75% of them evolved from progenitor ZAMS binaries with F-type or 
later primaries. This places the majority of the ZAMS binary progenitors of ZACVs with 
BP secondaries within the brown dwarf desert. Therefore, it is also possible that the lack of 
observed CVs with orbital periods shorter than the period minimum may be due to a lack of 
progenitor binaries from which to draw. Interestingly, if hydrodynamical models of the CE 
phase eventually do not indicate the existence of a cutoff mass > 0.1 M Q , this could provide 
strong, independent support for the existence of the brown dwarf desert. 
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4-2.2. The lower edge of the period gap: an imprint of CV formation''. 



Willems et al.l (120051 ) have computed comprehensive model populations of present-day 
CVs that formed with orbital periods less than 2.75 hr, including models with an additional 
source of angular momentum loss besides gravitational radiation. As an example of such a 
source, they considered circumbinary disks. Their model populations that include circumbi- 
nary disks successfully resolve the well-known dis crepancy betwe en t he predicted and the 
observed value of the period minimum in CVs (see |Pattersonlll998l and lKolb fc Baraffdll999 
for a discussion of this discrepancy). However, these models fail to reproduce the sharp lower 
edge of the period gap at 2.25 hr seen in the observed distribution. Willems et al. note that 
agreement could be improved by "a significant modification of the ZACV population near 
the lower edge of the period gap" or by a flow of systems into the gap from CVs formed 
above the gap. They simulate the latter flow by artifically multiplying the birthrate at 2.25 
hr in their models by a factor of ~ 100. However, as shown in a very recent follow-up paper 
(jWillems et al.l 120071 ) , this flow is i n consis tent with circumbinary disk models of CVs that 
form above the gap. I Willems et al.l (120071 ) find that the additional angular momentum loss 
created by the circumbinary disk increases the mean mass transfer rate in CVs formed above 
the gap and causes these systems to bounce at an orbital period of 2.75 hr (the upper edge 
of the gap). Therefore, CVs that formed above the gap in their models with circumbinary 
disks will always remain above the gap a nd, thus, canno t creat e the flow of systems to orbital 
periods shorter than 2.75 hr required by lWillems et al.l (120051 ). Consequently, it would seem 
that better agreement between their models and observations concerning the lower edge of 
the period gap must be achieved by a modification of the standard model for the ZACV 
population. 

In our PL model sequence for ZACVs, we find that as the value of n is increased, the 
fraction of CVs formed below the period gap decreases and the fraction of CVs formed in 
the gap increases (see Table 3). This trend is especially evident in the PL 2.0 model, where 
the fraction of CVs that form below the gap is reduced by over a factor of 7 and the fraction 
of CVs that form in the gap is increased by over a factor of 4 compared with our standard 
model. In addition, we find that as n is increased, the short-period peak in the orbital period 
distribution of ZACVs shifts to longer orbital periods (see Fig. 8). For our standard model, 
this peak occurs at an orbital period of 1.8 hr (log P = -1.125 in Fig. 6). For the PL 0.5, 1.0 
and 2.0 models, this peak occurs at an orbital period of 1.9 hr (log P = -1.10), 2.1 hr (log P = 
-1.05), and 2.4 hr (log P = -1.0), respectively. Interesti ngly then, this shift pr ovides precisely 
the modification of the ZACV population required by I Willems et al.l (120051 ). We note that 
in the CUTM model sequence, increasing the cutoff mass also shifts the short-period peak 
in the orbital period distribution to longer orbital periods, although to a lesser degree. For 
the CUTM 0.15 model, this peak occurs at 2.0 hr (log P = -1.175 in Fig. 9). Extrapolating 
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from our CUTM model sequence, for this peak to occur at 2.25 hr, a cutoff mass of 0.27 M Q 
would be required. 

We suggest, therefore, that it is possible that the sharp lower edge of the period gap 
may not be due solely to secular evolution, but could als o be, in par t, an imprint of CV for- 
mation. In the standard model of CV secular evolution (jKingill988l ). systems formed above 
the gap evolve to shorter orbital periods on a rapid timescale due to magnetic braking. 
Such rapid evolution eff ectively washes out any feature s from the underlying ZACV orbital 
period distribution (see iKolbl Il993l ; iHowell et al.ll200ll ). making such an imprint unlikely. 
However, recent observations of magnetic activity in single stars in young open clusters 



(ISills et al 



2000 



Pinsonneault. Andronov &: SiUsll2002l ) and of X-ray emission from rotating 



dwarfs ( Pizzolato et al. 20031 ) have called into question the m agnetic braking prescription 



used in the standard model (disrupted magnetic bra king, e.g.. iRappaport. Verbunt fc Joss 
1983 ). Magnetic braking prescri ptions developed by Andronov et al. ( 20031 ) based on the 
open cluster observations and by llvanova fc Taaml (120031 ) based in part on the X-ray emis- 
sion data produce angular momentum loss rates above the gap that are significantly lower 
than in standard magnetic braking. It is possible, therefore, that in a model of secular evo- 
lution where the rapid flow of systems from above to below the gap is reduced (as in the 
prescriptions by Andronov et al. and Ivanova & Taam) or eliminated (as in secular evolution 
with circumbinary disks), features in the ZACV orbital period distribution may persist and 
be evident in the present-day (evolved) CV distribution^ In order for the lower edge of the 
period gap to be such an imprint, the standard model of CV formation must be modified to 
produce a peak in the ZACV orbital period distribution at ~ 2.25 hr. One way of making 
such a modification is to assume that acE is a moderate to strong function of secondary mass 
(e.g., n > 1 in PL model sequence). However, we note that our preliminary investigation of 
the dependence of acE on the evolutionary state of the giant indicates that such a peak can 
be achieved in other ways as well. 



4.3. The dependence of acE on other orbital parameters 

It is very unlikely that cxce depends solely on secondary mass. As a simple example, if 
the efficiency of energy transfer depends upon the magnitude of the drag torque, then acE 
would not only depend upon the secondary mass, but also upon the density profile within 
the giant. This profile is determined by the core mass and total mass of the giant. Thus, 



6 We note that 
of the period gap. 



Webbink 



(|1979l ) made a somewhat similar argument in his so-called "static" explanation 
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atCE would be a multi-dimensional function of M s , M c , and M p in this scenario. In general, 
we would expect that a C£ is at least a function of secondary mass and the evolutionary 
state of the giant star, which is primarily determined by the giant's core mass. For example, 
acE might be different for CE evolution involvi ng a first giant branch primary than for CE 
evolution involving an AGB primary (see e.g., ISandquist et al.ll2000l ; IWillems et al.l 120051 ). 
To address the potential multi-dimensional nature of acE, we have chosen to investigate 
the dependence of oice upon secondary mass in this first paper and its dependence upon 
the evolutionary state of the giant in paper II. Clearly, we can only hope to approximate 
the true dependence of acE on the orbital parameters of the binary by considering its 
dependence on each parameter separately. Therefore, we urge caution when interpreting our 
results, particularly when it comes to inferring a particular q:ce{M s ) dependence based on 
comparison with observations of PCEBs and ZACVs. Observed trends in a given orbital 
parameter distribution for either population could also be due to strong dependences of acE 
on orbital parameters other than secondary mass, for example. 



4.4. Guiding detailed theoretical models of the CE phase and future 

observations of PCEBs 

One of our main goals in conducting the present study is to identify interesting regions of 
orbital parameter space that deserve further study by researchers who are calculating detailed 
hydrodynamical models of the CE phase. As mentioned in section 1, such calculations are 
extremely time intensive and demand an efficient use of available computational resources. 
To that end, we believe that calculating hydrodynamical models for a fixed, typical giant 
mass and a sequence of low mass secondaries should be given some priority. Our model 
populations of ZACVs indicate that if acE is very low (< 0.20), the distribution of ZACVs 
with orbital periods less than ~ 3 hr will be significantly altered compared with our standard 
model in which o>ce = 1-0 (see Fig. 6). Such a modification of the ZACV population may 
have important ramifications for our understanding of the orbital period distribution in 
CVs below and in the period gap. As mentioned in section 1, oice may be very small 
for CE evolution involving a very low mass secondary. In such cases, the rate of orbital 
energy dissipation within the CE is becoming comparable to the luminosity of the giant and 
insufficient to drive the bulk hydrodynamical motion needed the eject the envelope (R. Taam, 
2003, private communication). However, there have only been two detailed hydrodynamical 
calculations of CE evolution involving a secondary with M s < 0.30 and both involv ed 



giants with core masses < 0.45 M Q (i.e. on the first giant branch) (ISandquist et al.ll2000i ) 



We suggest calculating a sequence of models with an AGB primary that has a 0.6 M Q 
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core and secondaries ranging in mass from 0.05 M to 0.30 M . The number of models 
calculated will depend upon the available computational resources, but should include a 
sufficient number of models to estimate a dependence of cxce on secondary mass. At a 
minimum, models with M s = 0.05, 0.10, 0.20 and 0.30 M Q should be calculated. If the 
evolution involving the 0.10, 0.20 and 0.30 differ significantly, then additional models with 
0.15 M Q or 0.25 M may need to be calculated. 

If the above sequence of models indicates the existence of a cutoff mass for the secondary 
below which merger always occurs, then we suggest that additional CE models with that 
cutoff mass be calculated for other giant masses. If the value of the cutoff mass is found not to 
be strongly dependent upon the mass of the giant star, then such a cutoff mass may provide a 
prediction of hydrodynamical models of CE evolution that is feasible to test observationally. 
The chances of observing a binary in the CE phase are miniscule since the phase is so short- 
lived, making direct tests of detailed model calculations impossible. Consequently, we must 
rely on observations of PCEBs to provide indirect tests of hydrodynamical models of the 
CE phase. The value of a cutoff mass predicted from these hydrodynamical models does 
not require population synthesis calculations of present-day PCEBs in order to test it. Such 
population synthesis calculations are inherently dependent upon uncertain quantities such as 
the prescription for angular momentum loss following the CE phase (i.e., magnetic braking) 
and the mass ratio distribution in ZAMS binaries. Further, such calculations require a 
prescription for the CE phase, not just the value of an extremum for one orbital parameter. 
For these reasons, a predicted value of a cutoff secondary mass, assuming one exists, is 
clearly advantageous and would provide a direct, basic test of hydrodynamical models of CE 
evolution. 

Observationally, providing a sample of PCEBs that is both well-defined and statistically 
complete, especially with regard to very low mass secondaries, should be given some priority. 
Regardless of whether or not a cutoff mass for the secondary exists, such a sample is crucially 
needed if we are to have any confidence in the results of hydrodynamical calculations of CE 
evolution. 



5. Conclusions 

We have undertaken a systematic study of the effects of a CE efficiency parameter that 
is a function of secondary mass on models of the present-day populations of PCEBs and 
ZACVs and their distributions over orbital period, orbital separation, secondary mass and 
WD mass. We investigated three possibilities: (1) constant acE models for six values of acE' 
1.0, 0.6, 0.3, 0.2, 0.1, and 0.05, with acE = 1 designated as our standard model, (2) a power 
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law dependence, olce = {M s ) n , with n = 0.5, 1.0 and 2.0, and (3) a dependence in which 
o>ce approaches 1 for large secondary masses and acE = below some assumed cutoff mass, 
a CE = 1 - M cut /M s , with M cut = 0.0375, 0.075 and 0.15 M Q . We denoted these constant, 
power law and cutoff mass model sequences by CON, PL, and CUTM, respectively. 

The two most important findings of this study are: 

(1) Our CON model sequence for ZACVs shows that for values of cxce ^ 0.2, the 
predicted population of present-day ZACVs differs significantly from our standard model. 
All prior p opulation synthesis studies of the formation of CVs only considered valu es of 
a CE > 0.3 JPolitanolll988l . ll996l : [de~Koollll992l ; [Howell et alihoOll : Iwillems et al.ll2005h and 
found that their model populations were not strongly dependent upon the value of «cb- Our 
results indicate that a much wider range of values for ccce, including very low values, must 
be considered in order for a dependence to be seen. 

(2) Our model sequences for ZACVs in which acE is a function of secondary mass show 
that if otcE decreases rapidly for small secondary masses, the orbital period distribution in 
ZACVs for orbital periods in and below the gap differs significantly from our standard model. 
These differences are most evident in the PL 2.0 and CUTM 0.15 models. For these models, 
the fraction of CVs forming with orbital periods below the gap is reduced significantly, the 
fraction forming in the gap is increased significantly, and both the short-period peak and 
the minimum period in the ZACV orbital period distribution shift to considerably longer 
orbital periods. These results indicate that if o>ce is strongly depe ndent upon M ? for small 
secon dary masses, the standard model for the formation of CVs ( Ide Kooll Il992l ; iPolitano 
19961 ) with orbital periods below ~ 3 hr will need to be modified. 



Other findings of this study are: 

(3) Approximately 25% of all present-day ZACVs in our standard model are predicted 
to form with orbital periods shorter than 77 min, the observed orbital period minimum in 
CVs. These ultrashort-periods occur only in CVs that form with very low mass (late M) or 
substellar secondaries. Even with subsequent secular evolution, this prediction is in sharp 
conflict with the observed percentage of CVs with P < 77 min, which is < 1%. In our PL 
2.0 and CUTM 0.15 models, no CVs are formed with M s < 0.10 M because of the steep 
decline in olce for very low secondary masses. Consequently, no CVs are formed with P < 
77 min in these models, in much better agreement with the observations. This suggests the 
possibility that the observed scarcity of CVs with P < 77 min may be evidence that a cutoff 
secondary mass for CE evolution exists and is > 0.10 M Q . However, at the moment this 
is not clear, since this scarcity of ultrashort-period CVs could also be due to other factors 
such as a lack of CV progenitors with BD secondaries. If hydrodynamical models of the CE 
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phase eventually do not indicate the existence of a cutoff mass > 0.1 M , this could provide 
strong, independent support for the existence of the brown dwarf desert. 

(4) It is possible that the sharp lower edge of the period gap at 2.25 hr could be, in 
part, an imprint of CV formation. However, two key modifications of the standard model of 
CV formation and evolution are required in order for this to be true. First, the rapid secular 
evolution of systems to shorter orbital periods from above the gap must be significantly 
reduced or eliminated. This would require a modifi c ation o f standard magn e tic br aking 



similar to thos e proposed r ecently by lAndronov et al.l (120031 ) . Ilvanova &: Taaml (120031 ) . and 



Willems et al.l (120051 . 120071 ). Second, the short-period peak in the ZACV orbital period 
distribution must be shifted from ~ 1.80 hr to ~ 2.25 hr. Such a shift can be produced 
if o>ce is moderately to strongly dependent upon M s {n ~> 1 in our PL model sequence). 
Preliminary investigation of the dependence of oice on the state of the giant star, which we 
will report upon in paper II, indicates that this shift can produced in other ways as well. 

(5) The fraction of systems that survive the CE phase, f surv , in our CON model sequence 
decreases approximately as («ce) 0,5 from 0.50 to 0.07 as cxce is decreased from 1.0 to 0.05 
(see Table 1). In our PL and CUTM model sequences, f surv remains < 1/2 for all models 
and varies rather modestly from 0.39 to 0.48. 

(6) Over a wide range of values for ccce in our CON sequence, our models predict that 
the majority of present-day PCEBs contain secondaries with masses > 0.375 M , orbital 
periods > 1 day and orbital separations > 0.025 AU, with most having periods of ~ 3 days 
and separations of ~ 0.05 AU. These models further predict that the present-day population 
of PCEBs should contain roughly an equal number of systems with He and CO WDs if acE 
is globally near unity, but should be clearly dominated by systems containing CO WDs if 
cxce is globally small (< 0.30). 

(7) In our PL and CUTM model sequences for present-day PCEBs, the only distribution 
that varies significantly is the secondary mass distribution. In the PL sequence, as n is 
increased from 0.5 to 2.0, the fraction of PCEBs with low-mass secondaries decreases by 
a factor of ~ 4 from 0.26 to 0.06. The decrease in the number of low-mass secondaries is 
particularly severe for the n = 2.0 model. In this case, there are no present-day PCEBs with 
secondary masses less than ~ 0.10 M . In the CUTM model sequence, significant changes 
only occur near the cutoff mass and the distributions are nearly identical for M s > 0.375 
M . 

(8) In our CON model sequence for present-day ZACVs, if oice ^ 0.30, the majority 
of the systems contain He WDs and have orbital periods below the period gap. On the 
other hand, if qce ^ 0.20, then the majority of the systems contain CO WDs and there is 
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a significant shift in the orbital period distribution towards periods in and above the gap. 
The majority of ZACVs contain low-mass secondaries regardless of the value of a C E- 

(9) In our PL model sequence of present-day ZACVs, the majority of the systems contain 
He WDs regardless of whether the dependence of a C E on M s is weak (n = 0.5) or strong 
(n = 2.0). Systems with low-mass secondaries dominate if the dependence of acE on M s 
is weak to moderate (n < 1), whereas systems with high-mass secondaries dominate if the 
dependence is strong. Finally, systems with orbital periods below the period gap are favored 
if the dependence of acE on M s is weak, while systems in and above the period gap are 
favored if the dependence is moderate to strong. 

(10) In our CUTM model sequence of present-day ZACVs, the majority of the systems 
contain low-mass secondaries and have orbital periods below the period gap if M cut < 0.075 
M Q , while the majority of the systems contain high- mass secondaries and have orbital periods 
above the gap if M cut = 0.15 M Q . Further, approximately the same number of systems 
contain He WDs and CO WDs if M cut < 0.075 M , while the majority of the systems 
contain CO WDs if M cut = 0.15 M . 

(11) One of our main goals in conducting the present investigation is to provide guidance 
to those performing detailed hydrodynamical calculations of CE evolution. To that end, we 
suggest that calculating detailed models for a fixed mass giant and a sequence of very low 
mass secondaries, ranging from 0.3 M down to 0.05 M , should be given some priority. 
The results of our study indicate that understanding CE evolution involving very low mass 
secondaries may have important ramifications for our understanding of the formation of 
CVs in and below the period gap. In addition, if hydrodynamical models of the CE phase 
determine that a cutoff secondary mass exists and is not strongly dependent upon the core 
mass of the giant star, then the value of such a cutoff mass could provide a prediction of 
these models that is independent of any population synthesis calculations and is feasible to 
test observationally. 

(12) We further suggest that observing a well-defined sample of PCEBs that is statisti- 
cally complete, particularly with regard to very low mass secondaries, should be given some 
priority. Such a sample will provide important observational tests of hydrodynamical models 
of the CE phase. 

Warm thanks to C. Carden, O. De Marco, D. Drapes, B. Gansicke, D. Hoard, K. 
Simkunas, R. Taam, and B. Willems for very useful discussions about this work. We also 
thank the anonymous referee for comments that improved the quality and the clarity of the 
paper. This work was funded in part by NSF grant AST-0328484 and by a grant from the 
Wisconsin Space Grant Consortium to Marquette University. 
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Table 1. Fraction of systems, f surv , that survive CE evolution 



Model 


fsurv 


CON 1.0 


0.50 


CON 0.6 


0.41 


CON 0.3 


0.28 


CON 0.2 


0.22 


CON 0.1 


0.13 


CON 0.05 


0.07 


PL 0.5 


0.46 


PL 1.0 


0.43 


PL 2.0 


0.39 


CUTM 0.0375 


0.48 


CUTM 0.075 


0.46 


CUTM 0.15 


0.43 
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Table 2. Selected results for model populations of PCEBs 



Model 




Fraction with 




N/N stand 


He WDs 


CO WDs 


M 3 < 0.375 M 


P < 1 day 


CON 1.0 


0.51 


0.49 


0.36 


0.23 


1.00 


CON 0.6 


0.43 


0.57 


0.35 


0.25 


0.74 


CON 0.3 


0.31 


0.69 


0.33 


0.27 


0.45 


CON 0.2 


0.23 


0.77 


0.31 


0.27 


0.32 


CON 0.1 


0.10 


0.90 


0.26 


0.24 


0.17 


CON 0.05 


0.02 


0.98 


0.20 


0.19 


0.09 


PL 0.5 


0.45 


0.55 


0.26 


0.19 


0.80 


PL 1.0 


0.40 


0.60 


0.18 


0.15 


0.66 


PL 2.0 


0.34 


0.66 


0.06 


0.08 


0.50 


CUTM 0.0375 


0.51 


0.49 


0.33 


0.22 


0.92 


CUTM 0.075 


0.47 


0.53 


0.29 


0.20 


0.85 


CUTM 0.15 


0.43 


0.57 


0.21 


0.15 


0.72 



Table 3. Selected results for model populations of ZACVs 



Model 








Fraction 


with 






N/N atand 




He WDs 


CO WDs 


M 3 < 0.375 M 


P < 77 min 


P < 2.25 hr 


2.25 hr < P < 2.75 hr 


P > 2.75 hr 




CON 1.0 


0.58 


0.42 


0.72 


0.24 


0.67 


0.04 


0.29 


1.00 


CON 0.6 


0.54 


0.46 


0.72 


0.23 


0.65 


0.05 


0.29 


0.96 


CON 0.3 


0.46 


0.54 


0.69 


0.17 


0.58 


0.09 


0.34 


0.83 


CON 0.2 


0.41 


0.59 


0.68 


0.15 


0.53 


0.11 


0.35 


0.72 


CON 0.1 


0.26 


0.74 


0.61 


0.10 


0.41 


0.15 


0.44 


0.48 


CON 0.05 


0.12 


0.88 


0.51 


0.05 


0.27 


0.16 


0.56 


0.28 


PL 0.5 


0.46 


0.54 


0.67 


0.15 


0.58 


0.07 


0.35 


0.84 


PL 1.0 


0.31 


0.69 


0.57 


0.05 


0.42 


0.12 


0.46 


0.63 


PL 2.0 


0.08 


0.92 


0.35 


0.00 


0.09 


0.18 


0.72 


0.37 


CUTM 0.0375 


0.53 


0.47 


0.68 


0.18 


0.62 


0.05 


0.34 


0.86 


CUTM 0.075 


0.46 


0.54 


0.63 


0.05 


0.54 


0.06 


0.39 


0.74 


CUTM 0.15 


0.28 


0.72 


0.48 


0.00 


0.31 


0.13 


0.55 


0.52 
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Fig. 1. — Comparison of assumed functional forms for olce with estimated values of olce 
from detailed hydrodynamical calculations (crosses). The solid lines show a C E = (M s ) n for 
n = 0.5 (top curve) and n = 2.0 (bottom curve). The dotted lines show a C E = 1 — M cut /M s 
for M cut = O.O375M (top curve) and M cut = O.15M (bottom curve). 
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Fig. 2. — Theoretical present-day distributions of the orbital periods (top panels) and the 
orbital separations (bottom panels) in PCEBs for o>ce = 1.0 (black), 0.6 (magenta), and 0.3 
(blue) (left panels) and for a C£ = 0.2 (red), 0.1 (green) and 0.05 (purple) (right panels). 
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Fig. 3. — Theoretical present-day distributions of the secondary masses (top panels) and the 
WD masses (bottom panels) in PCEBs for o>ce = 1.0 (black), 0.6 (magenta), and 0.3 (blue) 
(left panels) and for acE — 0.2 (red), 0.1 (green) and 0.05 (purple) (right panels). 
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Fig. 4. — Theoretical present-day distributions of the orbital periods (top left), the orbital 
separations (top right), the secondary masses (bottom left) and the WD masses (bottom 
right) in PCEBs for olce — {M s ) n , where n = 0.5 (blue), 1.0 (green) and 2.0 (red). 
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Fig. 5. — Theoretical present-day distributions of the orbital periods (top left), the orbital 
separations (top right), the secondary masses (bottom left) and the WD masses (bottom 
right) in PCEBs for a CE = 1 - M cut /M s , where M cut = O.O375M (red), O.O75M (green), 
and 0.15M© (blue). 
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Fig. 6. — Theoretical present-day distributions of the orbital periods (top panels) and the 
orbital separations (bottom panels) in ZACVs for o>ce = 1-0 (black), 0.6 (magenta), and 0.3 
(blue) (left panels) and for a C£ = 0.2 (red), 0.1 (green) and 0.05 (purple) (right panels). 
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Fig. 7. — Theoretical present-day distributions of the secondary masses (top panels) and the 
WD masses (bottom panels) in ZACVs for o>ce = 1-0 (black), 0.6 (magenta), and 0.3 (blue) 
(left panels) and for acE — 0.2 (red), 0.1 (green) and 0.05 (purple) (right panels). 
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Fig. 8. — Theoretical present-day distributions of the orbital periods (top left), the orbital 
separations (top right), the secondary masses (bottom left) and the WD masses (bottom 
right) in ZACVs for olce — (M s ) n , where n = 0.5 (blue), 1.0 (green) and 2.0 (red). 
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Fig. 9. — Theoretical present-day distributions of the orbital periods (top left), the orbital 
separations (top right), the secondary masses (bottom left) and the WD masses (bottom 
right) in ZACVs for a CE = 1 - M cut /M s , where M cut = O.O375M (red), O.O75M (green), 
and 0.15M© (blue). 
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Fig. 10. — Observed orbital period distribution in CVs. The data were taken from the latest 
online version (RKcat7.6, 1 January 2006) of the iRitter fc Kolbl (120031 ) catalog. We do not 
include in the distribution any systems in the catalog that may be possible AM CVn systems. 



